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Abstract 


We present multiwavelength photometry and spectroscopy of SN 2022jli, an unprecedented Type Ic supernova 
discovered in the galaxy NGC 157 at a distance of ~ 23 Mpc. The multiband light curves reveal many remarkable 
characteristics. Peaking at a magnitude of g= 15.11 + 0.02, the high-cadence photometry reveals periodic 
undulations of 12.5 + 0.2 days superimposed on the 200-day supernova decline. This periodicity is observed in the 
light curves from nine separate filter and instrument configurations with peak-to-peak amplitudes of ~ 0.1 mag. 
This is the first time that repeated periodic oscillations, over many cycles, have been detected in a supernova light 
curve. SN 2022jli also displays an extreme early excess that fades over ~25 days, followed by a rise to a peak 
luminosity of Lop = = 10?*! erg s^! Although the exact explosion epoch is not constrained by data, the time from 
explosion to maximum light is Z 59 days. The luminosity can be explained by a large ejecta mass (Mej + 12 + 
6 Mc) powered by ? Ni, but we find it difficult to quantitatively model the early excess with circumstellar 
interaction and cooling. Collision between the supernova ejecta and a binary companion is a possible source of this 
emission. We discuss the origin of the periodic variability in the light curve, including interaction of the SN ejecta 
with nested shells of circumstellar matter and neutron stars colliding with binary companions. 


Original content from this work may be used under the terms 

BY of the Creative Commons Attribution 4.0 licence. Any further 
distribution of this work must maintain attribution to the author(s) and the title 
of the work, journal citation and DOI. 
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Supporting material: data behind figure, tar.gz file 


1. Introduction 


Stars with zero-age main-sequence masses (Mzams) greater 
than 8 Mo end their lives as core-collapse supernovae (CCSNe; 
Smartt 2009; Langer 2012), producing a diverse range of 
transients (e.g., Gal-Yam 2017; Modjaz et al. 2019). The variety 
in the observable properties of these SNe is thought to be 
dependent on the initial mass, metallicity, binarity, and mass- 
loss history of the progenitor star. Hydrogen-poor CCSNe are 
referred to as stripped-envelope SNe (SESNe) owing to 
significant mass loss of the progenitor, removing hydrogen, 
and in some cases helium, from the stellar envelope. SESNe 
classified as Type Ic do not show hydrogen or helium in their 
optical spectra, although the extent of helium stripping is still 
uncertain (Hachinger et al. 2012; Williamson et al. 2021). 
Envelope stripping can occur through strong stellar line-driven 
winds (e.g., Vink & de Koter 2005; Shenar et al. 2020) or 
interaction with a binary companion (Podsiadlowski et al. 1992). 

Evidence for periodicity has been searched for in SN light 
curves. Nicholl et al. (2016) investigated the undulations in the 
superluminous SN 2015bn but were limited by the duration of 
their time series and could not reliably identify periodicity. 
West et al. (2023) suggested a repeating pattern of 32 + 6 days 
in the declining light curve of SN 2020glb (an explosion 
somewhat similar to SN 2015bn). However, insufficient cycles 
were observed to perform robust statistical checks for 
periodicity. Martin et al. (2015) and Fraser et al. (2013) 
suggested a periodicity in the optical light curve of SN 2009ip, 
but Fraser et al. (2015) subsequently found no evidence for the 
periodicity in extensive R-band data. The light curve of the 
luminous, fast optical transient AT 2018cow was subject to 
periodicity searches, and while none were found in the optical, 
marginal evidence for periodicity in the variable X-ray light 
curve was suggested (Rivera Sandoval et al. 2018; Kuin et al. 
2019; Margutti et al. 2019). Perhaps the most promising 
detection of periodicity in SN emission is in the radio light 
curves of SN 1979C (Weiler et al. 1992) and SN 20011g (Ryder 
et al. 2004), which have been attributed to fluctuations in the 
density of the circumstellar medium (CSM) produced by binary 
stellar wind interactions. 

In this paper we present an extensive follow-up campaign of 
the Type Ic SN 2022jli from ~—50 to +200 days relative to 
maximum light. SN 2022jli presents an unusually long-lived, 
luminous early excess followed by a long rise time and slow 
spectroscopic evolution. The extensive, almost daily, photo- 
metric coverage of this bright SN for 200 days after peak 
indicates a periodic variability (P = 12.5 + 0.2 days) observed 
in multiple bands and instruments, with amplitude of order 196 
of the peak bolometric luminosity of the SN. 


2. Discovery and Classification 


Libert Monard discovered a transient in NGC 157 from 
Kleinkaroo Observatory and submitted the discovery report on the 
Transient Name Server (TNS) as AT 2022jli on 2022 May 5.17 
UT at an unfiltered magnitude ~14 mag (Monard 2022). With the 
ATLAS survey (Tonry et al. 2018b; Smith et al. 2020), we 
independently detected the object (internal name ATLAS22oat) 
on 2022 May 1641 (at o= 14.3 mag). The original TNS 


Discovery Report of Monard (2022) registered the object with 
an astrometric error of 14" 4 and our ATLAS transient server 
(Smith et al. 2020), which dynamically links to TNS 
discoveries, did not associate the two sources. The ATLAS- 
automated TNS registration triggered a new source (a= 
8769038 6 = —8°38668), a discovery report, and name (AT 
2022jzy). To prevent confusion, AT 2022jzy was manually 
removed entirely from the TNS records, and the original 
incorrect coordinates of AT 2022jli were replaced with those 
from ATLAS while preserving L. Monard's discovery credit 
(O. Yaron, private communication). A low-resolution (R= 
100) spectrum from a 0.35m telescope (Grzegorzek 2022) 
indicated a likely Type Ic. This classification was confirmed 
(Cosentino et al. 2022) by the extended Public ESO Spectro- 
scopic Survey of Transient Objects (ePESSTO+; Smartt 
et al. 2015). 

The SN is offset by 35"2 N and 1588 W from the center of 
its host galaxy NGC 157, which has a redshift of z = 0.0055. 
The kinematic distance on the NASA/IPAC Extragalactic 
Database (NED), from the recessional velocity (corrected for 
Virgo infall and assuming Hy=73kms 'Mpc_'), is 
D=23 +2 Mpc. Distance estimates from the Tully—Fisher 
and Sosies methods (e.g., Terry et al. 2002; Tully et al. 2013) 
have a large range from 11 to 29Mpc, and we adopt the 
kinematic distance D = 23 + 2 Mpc throughout the rest of this 
paper. The foreground Milky Way reddening is Ay — 0.1186 
(Schlafly & Finkbeiner 2011), and given the position of SN 
2022jli, some internal host extinction is likely present. We do not 
account for possible host extinction, but this is likely to be low 
owing to the lack of narrow NaI D absorption in the spectra 
(Poznanski et al. 2012), and our main results are not sensitive to 
this choice. The spectroscopy and photometry presented in this 
paper have been corrected for foreground galactic extinction, and 
the spectra have been shifted into the rest frame. We note that 
NGC 157 also hosted the Type Ic SN 2009em (Monard 2009). 


3. Observations 
3.1. Imaging and Photometry 


The first observations of SN 2022jli were reported to the 
TNS by L. Monard (Monard 2022). To our knowledge there 
are no pre-explosion nondetections available, as the object had 
just emerged from solar conjunction. Monard reported four 
epochs of unfiltered CCD photometry from observations taken 
at the Kleinkaroo Observatory between 2022 May 5 (MJD 
59704) and 2022 May 22 UT (MJD 59721). 

The Asteroid Terrestrial-impact Last Alert System (ATLAS; 
Tonry et al. 2018b) began observing at the position of SN 
2022jli on 2022 May 16 (MJD 59715) in normal survey 
operations. ATLAS is a quadruple 0.5 m telescope system using 
broad orange (o, 5600-8200 A) and cyan (c, 4200-6500 A) 
filters. The combined four-telescope system surveys the 
observable sky to a typical 5e depth of ~19 mag and a 
cadence of 1-2 days. ATLAS photometry and astrometry are 
calibrated with the all-sky reference catalog sky (refcat2; 


34 L, Monard corrected this a day later in the TNS Comment for AT 2022jli, 
but the TNS database coordinates remained in error. 


THE ASTROPHYSICAL JOURNAL LETTERS, 956:L31 (13pp), 2023 October 10 


Tonry et al. 2018a). ATLAS photometry for SN 2022jli was 
obtained by forcing photometry at the location using the 
ATLAS forced photometry server (Shingles et al. 2021) and 
adopting a 3c clipped nightly mean. 

The Zwicky Transient Facility (ZTF; Bellm et al. 2019) 
observed the field beginning on 2022 July 03 UT (MJD 59763), 
giving the object the internal name ZTF22aapubuy. ZTF 
photometry in both g and r bands was obtained from the ZTF 
public stream using the Lasair’ broker (Smith et al. 2019). 

Photometry from ASAS-SN (Shappee et al. 2014) beginning 
on 2022 May 09 UT (MJD 59708) was obtained using the 
ASAS-SN Sky Patrol website*® (Kochanek et al. 2017). We 
adopt a cutoff MJD of 59762, after which we do not include 
ASAS-SN g-band photometry in favor of higher signal-to-noise 
ratio ZTF g band. 

We triggered follow-up photometric observations of SN 
2022jli using the IO:O camera at the 2 m Liverpool Telescope 
(LT; Steele et al. 2004). Using the LT, we obtained six epochs 
of ugriz-band observations and an additional griz-band 
observation between MJD 59817 and 59894. 

The griBV-band photometry was obtained through the 
Global Supernova Project using the 1 m Las Cumbres 
Observatory (LCO; Brown et al. 2013). Additional V-band 
Observations were recovered from acquisition images taken 
with the ESO Faint Object Spectrograph and Camera version 2 
(EFOSC2; Snodgrass et al. 2008) on the ESO 3.58 m New 
Technology Telescope (NTT; Wilson 1983) during spectro- 
scopic follow-up by PESSTO (Smartt et al. 2015). 

All CCD reductions were performed using instrument 
specific pipelines. Photometric measurements for the LCO 1 
m, EFOSC2, and griz LT-IO:O data were made using 
AUTOPHOT (Brennan & Fraser 2022) without host subtraction. 
Photometry in griz bands was calibrated against Pan-STARRS 
field stars (Flewelling et al. 2020), and BV-band photometry 
was calibrated using the APASS catalog (Henden et al. 2016). 
LT u-band measurements were performed using the PSF 
package?” (Nicholl et al. 2023) and calibrated against the Sloan 
Digital Sky Survey (SDSS; Alam et al. 2015) catalog. 

The Gaia satellite (Gaia Collaboration et al. 2016), operated 
by the European Space Agency (ESA), observed SN 2022jli 
(internal name Gaia22cbu) between 2022 May 11 (MJD 
59710) and 2022 June 30 UT (MJD 59760). The Gaia Science 
Alerts Project (Hodgkin et al. 2021) reported three epochs of 
G-band photometry.’ We assume a pessimistic Gaia 
photometric uncertainty of 0.1 mag. 

Ultraviolet (UV) and optical photometry of SN 2022jli was 
performed with the Ultra-Violet and Optical Telescope 
(UVOT; Roming et al. 2005) on board the Neil Gehrels Swift 
Observatory (Swift, Gehrels et al. 2004). Swift observed the 
field 13 times between 2022 August 17 (MJD 59808) and 2022 
December 27 (MJD 59940) in the U, B, V, UVW1, UVM2, and 
UVW2 bands. The images in each filter were co-added, and SN 
magnitudes were extracted using standard tasks within the 
HEASOFT"? package. A small aperture of 3" was chosen, and 
an aperture correction was applied, following Brown et al. 
(2009). Without template subtraction, most UVW1, UVM2, and 
UVW2 exposures were nondetections. Keeping only detections 


35 https: / /lasair-ztf.lsst.ac.uk/object/ZTF22aapubuy / 
36 
https: / /asas-sn.osu.edu 
37 http: / [github.com/mnicholl/photometry-sans-frustration 
38 http:/ /gsaweb.ast.cam.ac.uk /alerts /alert/Gaia22cbu / 
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greater than the limiting magnitude, we retain only one epoch 
of UVM2 photometry but retain most observations in the UBV 
bands. The extinction-corrected light curve of SN 2022jli is 
shown in Figure 1. 


3.2. Radio Observation 


We obtained a single radio observation at the position of SN 
2022jli on 2022 September 19, starting at 22:15 UT with the 
enhanced Multi-Element Remotely Linked Interferometer 
Network (e-Merlin; DD14001; PI: Rhodes). Observations were 
obtained at a central frequency of 5.08 GHz with a bandwidth 
of 512 MHz. The observation consisted of 6-minute scans of 
the target interleaved with 2-minute scans on the phase 
calibrator (J0039—0942). The observation ended with a scan 
of the flux calibrator (J1331+3030) and the bandpass calibrator 
(J1407+2827). The data were processed using the e- MERLIN 
custom CASA-based pipeline (Version 5.8, Moldon 2021). 

The pipeline averages the data in both time and frequency 
space, flags the data for radio frequency interference, performs 
bandpass and complex gain calibration, and splits out the 
calibrated target field. We performed some further flagging and 
imaged the data within CASA. During the observation, two of 
the six antennas dropped out, which impacted the quality of the 
final image. We did not detect any radio emission at the position 
of SN 2022jli, and we measure a final rms noise of about 
52 uJy beam ' (and a 3c upper limit of 156 Jy beam). 


3.3. Spectroscopy 


We present our spectra spanning three epochs, —39 days to 
+47 days with respect to g-band maximum, and also show the 
low-resolution spectrum of Grzegorzek (2022) from the TNS 
(additional spectra will be presented in a separate publication). 
Foreground Galactic reddening was corrected using the 
dust extinction package of Astropy following the 
Fitzpatrick (1999) reddening law. All spectra presented in this 
work will be made available via the WISeREP repository 
(Yaron & Gal-Yam 2012). 

Follow-up spectroscopy was acquired using the EFOSC2 at 
the 3.58m NTT (Snodgrass et al. 2008) at two epochs through 
ePESSTO-+. The first EFOSC2 spectrum was taken on 2022 
May 24.42 UT (MJD — 59723.42), and our final spectrum was 
taken on MJD — 59811.15. The EFOSC2 grism used for the 
spectral sequence was Gr#13 (3685-9315 A). Data reductions 
were performed using the PESSTO pipeline, which includes 
flat-fielding, bias-subtraction, wavelength and telluric correc- 
tion, and flux calibration as described by Smartt et al. (2015). 

One epoch from the University of Hawaii 2.2 m telescope 
was obtained on 2022 July 24.62 (MJD 59784.62) using 
SNIFS (Lantz et al. 2004). The SNIFS spectrum was reduced 
using the Spectroscopic Classification of Astronomical 
Transients (SCAT) survey pipeline (Tucker et al. 2022). 

The extinction-corrected spectra are shown in Figure 2. 


4. Analysis 
4.1. Light Curve 


SN 2022jl fades for ~25 days after discovery in the go 
bands; using a linear fit, we measure a decline rate of ~5 mag 
(100 days) !, which is incompatible with ?Co decay. This 
early decline is not well sampled with low-cadence coverage 
from only ASAS-SN (g), ATLAS (o), Gaia, and a single V- 
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Figure 1. Top: multicolor extinction-corrected light curves of SN 2022jli including 


photometric errors. Black lines along the bottom of the plot indicate spectroscopic 


observations. Open symbols are synthetic photometry performed on the EFOSC2 spectrum. Bottom: light curves after detrending with a second-degree polynomial fit 
including only the region after maximum light (vertical dashed line). Overplotted for each band are the best-fit sinusoids produced in the GLS analysis (Section 4.2). 


We retain the colors from the top panel and remove variation of point shape and do 
MJD 59750. 


(The data used to create this figure are available.) 


band observation from the NTT. The light curve begins to rise 
after MJD 59734, reaching maximum light on MJD 59763, 
indicating that the main peak is at ~59 rest-frame days from 
discovery. SN 2022jli exhibits a significantly longer rise to 
maximum light than literature samples of SESNe (Prentice 
et al. 2019), comparable to a small subset of slowly evolving 


not include photometric errors for visual clarity, including only observations after 


Type Ibc SNe (e.g., Anupama et al. 2005; Lyman et al. 2016; 
Taddia et al. 2016, 2018; Karamehmetoglu et al. 2023). The 
SN eventually fades in the optical at a rate of —1 mag 
(100 days) !, which is compatible with with *°Co decay, 
suggesting a radioactively powered main peak (Woosley et al. 
1989). The r band declines 0.4 mag (100 days) ' faster than the g 
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Figure 2. Spectroscopic follow-up observation of SN 2022jli. Spectra are labeled by time since maximum light in the rest frame. Common SN lines are marked as 
vertical lines corresponding to the rest wavelengths. We include here the Grzegorzek (2022) classification spectrum. A TARDIS (Kerzendorf & Sim 2014) model of 
H- and He-free material is also included. Spectra for the Type Ic SNe SN 2004gk (Shivvers et al. 2019) and SN 19941 (Modjaz et al. 2014) were obtained from 


WISeREP (Yaron & Gal-Yam 2012). 


band, showing an unusual evolution toward bluer colors in g — r 
over time. The double-peaked early light curve indicates that Ni 
decay cannot account for the full structure of the light curve. 

As the light curve slowly fades from peak, the extensive high- 
cadence photometry of SN 2022jli captures clear undulations in 
the photometry. The undulations are visible in multiple bands (B 
to 7) and across different telescope and instrument combinations, 
indicating that this is neither an instrumental nor a calibration 
effect. We subtract the SN continuum and reveal these 
undulations more clearly in the bottom panel of Figure 1. We 
identify repeating bumps with a consistent timescale for all 
filters and discuss this in detail in Section 4.2. 

During the initial decline, we observe an increase in the 
ASAS-SN g-band photometry on MJD 59722.37. Seeking to 


confirm the validity of this observation, we perform synthetic 
photometry on the EFOSC2 spectrum taken on MJD 59723.42, 
which was calibrated to the V-band acquisition image. We 
include the synthetic photometry in Figure 1; the errors on the 
points are consistent with the ASAS-SN g-band photometry 
and other contemporary photometric observations in G and o 
bands. We interpret this epoch as a short-lived, luminous 
episode during the initial excess. 


4.1.1. Bolometric Light Curve 


We compute a pseudobolometric light curve by integrating 
under the BgcVrGoiz-band observations using the publicly 
available code SUPERBOL (Nicholl 2018). From our 
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bolometric light curve we measure a peak luminosity 
Lop = 1072085004 erg s~!, which is within the typical range 
of SESNe found by Prentice et al. (2019). The total integrated 
luminosity (across the wavelength range covered by our filters) 
is Eop © 2.5 X 10” ergs. 

We compare the (pseudo)bolometric light curve to other 
SNe, including normal SESNe and those with double-peaked 
light curves, in Figure 3. SN 2022jli exceeds the peak 
brightness of the Type Ic SN 2007gr (Hunter et al. 2009), the 
Type Ib SN 2008D (Soderberg et al. 2008; Modjaz et al. 2009), 
and the relatively faint Type Ib SN 2007Y (Stritzinger et al. 
2009), showing a significantly more luminous broad peak and 
slower decay. SN 2007gr and SN 2007Y both display a 
monotonic rise and smooth decline, typical of a normal Type 
Ibc SN, unlike the double-peaked structured light curve of SN 
2022jli. The overall shape of the light curve resembles the 
unusual Type Ic iPTF15dtg (Taddia et al. 2016). Both SNe 
have a fast-declining early excess with a broad persistent 
maximum. SN 2005bf (Anupama et al. 2005) has an early peak 
and broad maximum but is significantly more luminous than 
SN 2022jli and declines significantly faster. 


4.1.2. Light-curve Modeling 


We model the light curve using simple models to derive a 
representative ejecta mass estimate for SN 2022jli using the 
Modular Open Source Fitter for Transients (MOSHIT; 
Guillochon et al. 2018). MOSHIT is a publicly available code 
that we use to fit semianalytic models to the multiband 
observed light curves of SN 2022jli. We use two models, one 
where we model only the broad “main” peak assuming 
radioactive decay of ??Ni as the only energy source (Arnett 
1982; Nadyozhin 1994), and another where we fit the full light 
curve interpreting the initial excess as shock cooling emission 
(SCE) from interaction with a CSM and a subsequent 
radioactively powered “main” peak (Chatzopoulos et al. 
2013). All model fitting was performed using the dynamic 
nested sampler DYNESTY package (Speagle 2020) option in 
MOSHT. 

The bottom panel of Figure 3 shows the ?9Ni-only model fit 
to SN 2022jli. To construct this model, we used the nickel- 
driven explosion model built into MOSFiT (Nadyozhin 1994), 
omitting any data during the early excess (before MJD 59732). 
We modify the priors of the model to require the explosion 
time to be before discovery, i.e., MJDexptosion < MJDaiscovery: 
The opacity was fixed at x = 0.1 cm? g and &,— 0.027 cm? 
g '. With no prior on ejecta velocity, the data require 
Vej & 2500 kms ! and M.j~4 Mo; this velocity is much 
lower than the measurements of the Fe II lines in the spectra 
(see Section 4.3). The posterior distribution for this fit is 
included in the Appendix (Figure 5). 

This model reproduces the maximum luminosity for the g 
and c bands but fits poorly to the redder bands, underestimating 
the flux particularly in the o and i bands. This simple model 
also fails to reproduce the fast g-band rise after the light-curve 
dip and the o-band peak luminosity. Color differences between 
the model light curve and the observed data are likely due to 
the blackbody assumption made by MOSFiT, as the true 
spectrum is dominated by strong emission and absorption lines 
by the time of maximum light. Further detailed modeling is 
warranted using more sophisticated techniques, which is 
beyond the scope of this work. For comparison, we perform 
an additional Arnett model (Arnett 1982) fit to the bolometric 
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Figure 3. Top: pseudobolometric light-curve comparison of SN 2022jli and 
other SNe with prominent bumps and well-studied Type Ibc SNe. The 
bolometric light curves of the Type Ic SN 2007gr (Hunter et al. 2009), Type Ib 
SN 2007Y (Stritzinger et al. 2009), Type Ic iPTF15dtg (Taddia et al. 2016), 
Type Ib SN 2008D (Soderberg et al. 2008; Modjaz et al. 2009; pseudobolo- 
metric light curve calculated in Nicholl et al. 2015), and Type Ib/c SN 2005bf 
(Anupama et al. 2005) were also constructed using SUPERBOL. The dashed 
line represents an Arnett model (Arnett 1982) fit to the SN 2022jli light curve 
after MJD 59732. Bottom: MOSFIT (Guillochon et al. 2018) multiband fit 
results for a °°Ni explosion model (solid lines). We include the best scoring 
realization from a CSM + *°Ni model fit to the first 110 days of SN 2022jli as a 
dashed line. 


light curve with a xX? fitting approach. Fixing & = 0.1 cm? g! 


and vej ~ 3000 kms !, we require Mii z 6 Mo, which is in 
agreement with the results from MOSFIT. This fit is included in 
Figure 3. 

To gauge the systematic modeling uncertainties within 
MOSFiT, we run the same nickel-driven model for the main 
peak with a range of vej to determine a probable range of Mej. A 
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fixed ejecta velocity of ve; = 3500 km s^! requires Ma £7 Mo, 
Vey = 6000 km s! requires Maz 18. Mo, and v4-— 
7000 kms! forces Ma47:21 Mc. With no pre-explosion 
nondetections available to constrain the explosion time and 
large systematic errors on the model, we adopt an indicative 
mass range for SN 2022jli of M, ~ 12 +6 Mo. 

In the second scenario we consider the contribution of SCE 
following shock breakout of the ejecta through a dense CSM 
using the *°Ni + CSM (CSMNI) model (Chatzopoulos 2013; 
Villar et al. 2017; Jiang et al. 2020). SCE is a natural 
interpretation for a fast-declining excess shortly after explosion 
and is now regularly detected for a range of SN subtypes. The 
same assumptions are made for the opacities as before, but the 
explosion time is left as a free parameter of the fit. We modify 
the code so that the interaction begins when the ejected material 
reaches the inner radius of the CSM at Ro. Our results indicate a 
CSM radius Ro~ 1 au, Ma ~ 20 Mo, Mcsm ~ 26 Mo, and 
Jni~ 0.01 for the early excess to be powered by CSM 
interaction. Our model (fit to the first 110 days) results in a 
poor match to the observed data, particularly the colors of the 
early excess, and requires physically improbable parameters; 
this model is shown in Figure 3. 

Although the Chatzopoulos et al. (2013) model implemen- 
ted in MOSFiT is relatively simple, the very large ejecta 
masses required imply that this scenario is physically unlikely. 
We return to this point in Section 5.2. We also emphasize 
that there is no robust measurement of explosion time 
to constrain the model since the earliest epoch from the 
Monard observations is after the SN appeared from solar 
conjunction. 

A long rise time (exceeding 30 days) is a rare occurrence for 
Type Ibc SNe (Lyman et al. 2016): long-duration light curves 
with rise times similar to SN 2022jli arise from only ~6%-10% 
of SNe Ibc in a bias-corrected sample (Karamehmetoglu et al. 
2023). We estimate that Maj = 122-6 Mo is required to 
provide the long rise to maximum light. An ejecta mass this 
extreme is rare and points to a high-mass progenitor star. 


4.2. Periodic Variability 


The declining light curve (Figure 1) shows ~0.05 mag 
undulations, which are present across all bands and appear to 
repeat with a regular amplitude and period. To search for and 
quantify any periodicity, we first removed the decline signature 
from the light curve. From the post-peak light-curve data 
(MJD > 59760), we produced a residual light curve in each 
band by fitting and subtracting a fourth-order polynomial fit 
(the lowest order that removes the SN decline) between MJD 
59760 and the end of the time series for each band. We applied 
this method to each of the BgcVroi bands independently and 
include the results in the bottom panel of Figure 1. The residual 
light curves show consistent oscillations over time across all 
bands. The periodicity in each of the residual light curves was 
quantified by computing a periodogram using a generalized 
Lomb-Scargle (GLS) method (Zechmeister & Kürster 2009). 
The periodograms for each band and a phase-folded light curve 
are shown in Figure 4. 

We find that the undulations have a dominant frequency of 
~0.08 day ! (or a period of ~12.5 days), where significant 
power is observed across the period range of 12-13 days, 
which is safely below the At/3 cutoff adopted by Martin et al. 
(2015) and Nicholl et al. (2016), where At is the length of the 
time series. The maximum GLS power in this region exceeds 
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Figure 4. Top: GLS periodograms from BgcVroi bands, where dashed 
horizontal lines correspond to the 196, 0.196, and 0.0196 FAP levels. The 
bottom panel shows a Lomb—Scargle periodogram generated from all bands 
simultaneously using gatspy. For the gr bands additional periodograms are 
computed using only LCO and ZTF photometry; the associated FAP levels for 
these periodograms are not plotted. All periodograms regardless of band or 
telescope—instrument combination show a significant power at ~12.5 days. 
Bottom: binned phase-folded ZTF g-band and ATLAS c-band light curve. 
Each band has been folded on their respective best periods; the unbinned 
phase-folded points are shown as open symbols. 


the 0.0196 false-alarm probability (FAP; see Zechmeister & 
Kürster 2009) level in all bands (shown in Figure 4). The 
detrended data reveal peaks and troughs with amplitudes in the 
range of 0.04-0.08 mag across all bands. Motivated by the 
synchronized behavior of the multiband photometry, we 
compute a periodogram fitting all the BgcVroi photometry 
simultaneously using the package gatspy (VanderPlas & 
Ivezié 2015; VanderPlas 2016), which generates a Lomb 
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—Scargle (Scargle 1982) periodogram for a multiband time 
series (Figure 4). The best period for the combined multiband 
observations is ~ 12.5 days. 

The consistency of the observed periodicity across the time 
series was verified through empirical mode decomposition (EMD; 
Huang et al. 1998, 1999), which is ideally suited to oscillatory 
detections in the presence of nonlinear and nonstationary 
processes that may impact the light curves of SNe. Following 
the methodology outlined by Jess et al. (2023), the intrinsic mode 
functions (IMFs) are extracted from the detrended time series. 
Subsequently, a Hilbert-Huang transformation (Huang & Wu 
2008) was performed to investigate the instantaneous frequencies 
across the observing window. A low-order IMF exhibits a 
frequency associated with a ~ 12.5-day period for the majority of 
the time series with little variation. Hence, the EMD processes 
applied here directly and independently support the GLS 
periodograms depicted in Figure 4. 

We recalculate the bolometric light curve using only the 
gcroiz bands to avoid washing out the periodicity with 
interpolation. Following identical detrending methods to the 
bolometric light curve, we measure the size of the oscillations 
(peak to trough) to be ~2 x 1079 erg s! over the underlying 
radioactively powered flux, which is on the order of 1% of the 
peak bolometric luminosity. We perform numerical integration 
under a single bump in the bolometric light curve to estimate 
the radiated energy Eyad bump ~ x 10% ergs. 


4.3. Spectra 


The spectroscopic evolution of SN 2022jli is shown in 
Figure 2, spanning from —53 days before maximum light to 
+48 days after. This includes a spectrum during the first 
maximum, which is unusual for SNe with a short-lived early 
peak. The line identifications in this section are based on those 
by Hunter et al. (2009). 

The first spectrum obtained during the early excess displays 
P Cygni absorption features of NaI D AA5891, 5897 and strong 
Fe ll AM924, 5018, 5169 absorption, typical of Type Ic SNe. 
The broader spectral coverage of the EFOSC2 spectrum (—39 
days) reveals Call H and K lines and the Ca II near-IR triplet. 

The spectra from +21.5 days show the emergence of a Sc II 
feature and have a complex blend of narrow emission lines and P 
Cygni features like some interacting SNe. The forbidden [Ca 1] 
AA7291, 7324 lines are prominent in the later spectra. We 
measure the velocities using Gaussian fits to the Fell \\4924, 
5018, and 5169 P Cygni absorption troughs as a proxy for the 
photospheric velocity. At —39 days we measure 8500+ 
300 km s !, 7000 2- 300 kms ' at +21 days, reducing further 
to 6700 + 300 km s^' by the +48-day spectrum, representing a 
slow recession of the photosphere inside the ejecta. 

We note the resemblance of the SN 2022jli spectra to SN 
1994] but also to SN 2004gk, although the spectroscopic 
evolution is not analogous to either object. Given the broad 
complex at Ha, which may be contaminated by Si II and C II, 
we cannot rule out H in the ejecta or CSM shell(s), but from the 
line ratios of other H lines we expect the contribution of H to be 
small. We also fail to identify unambiguous signatures of He in 
the spectra. 

We include comparison spectra of two representative Type 
Ic events, SN 2004gk (Shivvers et al. 2019) and SN 1994I 
(Modjaz et al. 2014). The post-peak (--25 days) spectrum of 
SN 2004gk best resembles the —39-day spectrum of SN 
2022jli. This suggests that SN 2022jli has undergone an 
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unusual evolution to resemble a post-peak normal Type Ic 
spectrum at this phase. 

We produced a TARDIS (Kerzendorf & Sim 2014; 
Kerzendorf et al. 2023) model (Figure 2) with the aim of 
reproducing the main spectral features of SN 2022jli during the 
first EFOSC2 spectrum. The model is a simple uniform 
abundance model with an H- and He-deficient composition 
dominated by C, O, Si, and Mg and a photospheric velocity of 
7500 kms !. We adopt a fexpiosion parameter (time since the 
start of homologous expansion) of 42 days before the 
observation at —39 days to best match the observed features. 
We successfully reproduce the prominent Fe II and Ca II features 
and continuum shape and show that a plausible Type Ic SN 
composition can reproduce the spectrum. The model does not 
reproduce the bump at 6500 A or the emission at 6150 A. The 
TARDIS configuration file is available in a .tar.gz package. 


5. Discussion 


The data presented in this paper show that SN 2022jli is 
unusual in many respects. The duration of the initial excess 
(225 days with no constraining nondetection) is unprecedented 
for a Type Ic SN. The bolometric light curve peaks at least 59 
days after explosion and could be longer given the uncertainty 
in explosion epoch. In combination with the periodic 
undulations, the SN 2022jli observational data set is unique. 


5.1. Scenarios for Periodic Variability 
5.1.1. Interaction with CSM 


The bumps we observe in the light curve could be due to 
ejecta interacting with concentric shells of circumstellar 
material. During the undulation, SN 2022jli is overluminous 
by ~1 x 10*° erg s^! for 12.5 days. Using the scaling relation 
L= 5Mcsuv? Ji rise (Smith & McCray 2007; Quimby et al. 
2007; Nicholl et al. 2016), we estimate the mass required for 
each bump to be Mcsm,bump © 10? Mo, assuming v= 
7000kms! from Fell line velocity measurements, and 
tise = 6.3 days. 

The average pre-explosion mass-loss rate needed to produce 
this CSM mass per undulation can be calculated by setting 
M [vw = Mcsm,bump/AR, where v, is the wind velocity and 
AR = Vfpump is the radial distance bounding this CSM mass. 
For an SN velocity v = 7000 kms ^ ! and foump = 12.5 days, this 
gives M = few x 1075(v,,/1000 km s7!) Mo yr ^l, This is 
consistent with the mass-loss rates observed from Wolf— Rayet 
stars (e.g., Sander & Vink 2020), suggesting that a "typical" 
wind mass-loss rate could potentially provide the CSM 
structure needed to explain the periodic undulations of SN 
2022jli, if subjected to a periodic modulation. 

Nested shells of dust caused by colliding winds in a massive 
binary system have recently been spectacularly revealed in 
JWST imaging (Lau et al. 2022). They showed that the 17 
observed shells were due to repeated dust formation episodes 
every 7.93 yr modulated by periastron passage of the companion 
O5.5fc star in the mutual orbit around the WC7 Wolf—Rayet 
star. An ejection velocity of v = 7000 km s^! means that the SN 
shock front travels Ar ~ 54 au in 12.5 days. By comparison, the 
nested dust shells around WR 140 are Ar = 4380 + 120 au. If 
SN 2022jli undulations were due to peaks in CSM density 
similar to WR 140 shells, a binary progenitor would need to 
eject these shells on timescales ~100 times more frequent, or 
with a ~0.2 yr periodicity. We note that the dust emission in the 
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shells in WR 140 does not necessarily imply enhanced gas 
densities (Pollock et al. 2021). 

In this scenario of concentric shells or rings, one might 
expect that light-travel time effects could broaden the 
undulation timescale as the shock expands. As the SN ejecta 
hits the back and front of the shells at the same time, the light- 
travel time from the back to front increases as At, ~ 2Vejtexp /c, 
or about 9 days (after 200 days of expansion). While there will 
be an integrated signal from all parts of the shell, the effect 
should be to broaden the timescale of the undulations; should 
the ejecta be photospheric, the broadening effect will be less. 
The data do not clearly support such broadening. 

We can explore potential progenitor candidates using the 
BPASSV2.2.2 predictions (Eldridge et al. 2017; Stanway & 
Eldridge 2018), restricting our search to Type Ic progenitors 
(see Stevance & Eldridge 2021): we select hydrogen-deficient 
systems with surface mass hydrogen «0.0001 Mo and surface 
hydrogen mass fraction «0.01. We also restrict our search to 
helium-depleted SN progenitors, with helium mass fraction 
«0.3. We then look for systems that satisfy the period estimate 
of a WR 140—like scenario by searching for P = 0.070.015 
yr (the error is chosen to give a window of roughly +/—5 
days). Finally, we also impose a luminosity and temperature 
constraint (log(L/Lo) > 5 and log(T/T>) > 4.3), as we are 
looking for WR4-O star systems. We find 27 BPASS models at 
solar metallicity that fulfill these requirements, and including 
the initial mass function weighting, we would expect about 15 
such systems to be formed per 1 million Mọ. All these systems 
have primary star (SN progenitor) masses in a rather narrow 
range of 10—13 Mo, while the secondaries range from 24.5 to 
60 Mo. Although all these predicted systems have stellar winds 
around 0.8 x 10? M. yr ', similar to the estimated 
requirement to create the CSM as mentioned in Section 5.1.1, 
a key factor in the periodic modulation is the eccentricity of the 
system. Stellar evolution models such as BPASS assume 
circularized orbits, so we cannot assess how many systems 
would be born with and maintain sufficient eccentricity. 


5.1.2. Accreting Compact Object 


An alternative mechanism to produce light-curve bumps in 
Type Ibc SNe was suggested by Hirai & Podsiadlowski (2022). 
After an SESN in a binary system, the newly born neutron star 
(NS) may receive a kick in the favorable direction of its 
companion. As a result, the NS may penetrate or skim the 
surface of the binary companion. They predict that material 
captured from the companion settles around the NS and that the 
accretion rate is likely to be super-Eddington. The accretion 
could result in outflows or jets that add further energy to the SN 
ejecta and result in additional luminosity. Accretion resulting in 
jets has been modeled by Hober et al. (2022), which they 
propose could power bumps in the late-time light curves of 
SNe. Undulations in the light curves of SLSNe have been 
detected (Nicholl et al. 2016; Inserra et al. 2017; Gomez et al. 
2021; Hosseinzadeh et al. 2022; West et al. 2023), but no 
repeating signature has been confirmed over multiple cycles. 

Hirai & Podsiadlowski (2022) calculate that even if only 
~0.01 Mois captured by the NS and if only ~1% of that is 
accreted, then the energy available would be of order 
Face © ace ov 10°° erg, which is comfortably enough to 
power a few percent of the total integrated SN flux of 
Eja 722.5 x 10°” erg. The direct interaction invoked in Hirai & 
Podsiadlowski (2022) requires a fairly fine-tuned kick direction 
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and velocity. However, a milder interaction as discussed in 
Hirai et al. (2018) and Ogata et al. (2021) may be sufficient. 
The companion star is inflated by heating from the SN ejecta 
—companion interaction, and the inflated part of the envelope 
may interact with the NS, causing periodic accretion on the 
timescale of the orbit. 

Accretion-powered jets after core collapse also have 
sufficient energy (Soker 2022; Hober et al. 2022) to power 
the excess flux observed in the light curve of SN 2022jli, but a 
modulation process is required. In the Hirai & Podsiadlowski 
(2022) scenario, the gradual inspiral of the NS into the 
companion is a pathway for the formation of a Thorne 
—Zytkow object (TZO). TZOs, which are NSs inside an 
envelope of nondegenerate diffuse material, have been 
predicted in the literature (Thorne & Zytkow 1975, 1977), 
but very few real candidates exist (O'Grady et al. 2020). An 
issue of the Hirai & Podsiadlowski (2022) scenario is that the 
orbit of the NS will decay rapidly (within ~5 orbits), making 
15 orbital cycles problematic. However, in the inflated 
companion case, the low-density envelope results in slower 
orbital decay (Hirai et al. 2018; Ogata et al. 2021). 

The accreting compact object model can be thought of as an 
internal powering source. The energy released must diffuse out 
through the ejecta on timescales determined by the opacity, 
density, and radius of the optically thick material. Hosseinzadeh 
et al. (2022) propose that a central origin is disfavored if the 


e : $ thump Af j 
dimensionless depth of the powering source, 6 c 2*7 9 is 


tise 
significantly less than unity. With Afjump = 12.5 and D ~ 60, 
this parameter then ranges between 0.2 and 0.7 for the earliest 
and latest bumps. This would marginally disfavor a central, 
internal powering source, although Hosseinzadeh et al. (2022) 
note that the expression is quite approximate and should only be 
treated as an order-of-magnitude result. 


5.2. Scenarios for Initial Maximum 
5.2.1. Interaction with CSM 


MOSFIT modeling of the early excess with an interaction- 
powered model requires Mj œ~ 12 Mo, which is compatible 
with the mass required for the long rise to maximum light. 
However, the duration requires significant Mcsw (^3 Mo in or 
modeling), which is very large for a Type Ic SN and would 
require an exotic mechanism to drive extreme mass loss shortly 
before explosion, such as pulsational pair-instability SN 
ejections (Woosley 2017). 

Perhaps the apparent duration of SCE is extended owing to 
enhanced opacity caused by Thomson scattering in the CSM 
(Moriya & Maeda 2012), which is eventually overtaken by the 
forward shock. Finally, the red spectrum during the first peak 
(Grzegorzek 2022) would appear inconsistent with luminous 
circumstellar interaction at early times (ie. with com- 
pact CSM). 

We cannot exclude CSM interaction as the source of the 
early excess; to do so would require more data to constrain the 
excess or more sophisticated modeling of the interaction to 
determine the viability of this scenario. 


5.2.2. Companion Collision 


Here we consider the emission from the collision of ejecta 
with the binary companion of SN 2022jli using the model 
suggested by Kasen (2010). In this scenario the interaction 
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shocks the SN ejecta, dissipating kinetic energy causing bright 
optical/UV emission. This additional contribution to the 
observed luminosity exceeds the radioactively powered SN 
for a short period, resulting in an early excess. We investigate 
the viability of this model using Equations (22) and (23) from 
(Kasen 2010) to estimate the luminosity and the collision 
luminosity timescale (f.) where (Leiso > Lyi) 


1/4 7/4 ,.—3/4 ,-1/2 E 
Leiso = WOM aia vdt dl a erg sl, 


(1) 


2/5 441/2,,3/10 ,.1/10 &Ni2/5, 4.2/5 
te < 713a M vo” «V us Mi dos 
"e 


(2) 


where a3 =a/ 10? cm (a is the orbital separation), M. — 
M/Ma is the ejecta mass M in units of the Chandrasekhar mass 
(Mcp), vo = v/10° ems "vw — 6 x 108 (E5,/M;)"? ems! 
(following Kasen 2010, we adopt ¢, = 1.69), Ky; is the opacity 
in the °°Ni-dominated region, and Ke is the ejecta opacity 
outside this region. The time since explosion is given as fgay, 
Mwios = Mwi/0.6 Mo, Es; = E/10?! ergs !, and E is the 
explosion energy. We adopt an indicative ejecta mass of 12 Mo 
from MOSHIT modeling in Section 4.1.2 and set ke = 
Kw; = 0.1 cm? g We set My; + 0.7 (for fni œ% 0.06 from 
MOSFIT) and set v, = 8500 km s ! from direct measurement 
of the —29-day spectrum (during the early excess). 

To produce the observed early luminosity on the order of 
^10? erg s ! and timescale of ~10 days, we would need to 
have separation ~1lau. For these parameters we calculate 
te S 17.5 days and Le iso 78 x 10% erg s^! for ty = 2 days and 
Loio£22 X 10 erg s ! for tj = 20 days. 

With these assumptions for separation and ejecta mass, our 
observations are compatible with the direct collision of the 
ejecta with the companion star. It is important to note that this 
scenario requires a favorable viewing angle; Kasen (2010) 
predicts that the collision should be visible in only ~10% of 
cases and that an orbital separation of ~1 au at the time of 
interaction may require the object to be close to pericenter. An 
important additional qualification of the Kasen (2010) model 
calculation is the assumption that the companion to the 
exploding star is filling its Roche lobe, as in a thermonuclear 
binary star explosion. Therefore, this calculation should be 
regarded only as an illustrative estimate of the energetics of a 
companion interaction. The progenitor systems considered in 
Section 5.1.1 have separations of 0.5 au and upward, which is 
compatible with the lau separation adopted for this 
calculation, although we note that they are typically not filling 
their Roche lobe and interacting. Should the excess be powered 
by companion collision, one might expect to observe late-time 
Ha emission from the companion and at late times observe a 
surviving but inflated companion star. 


6. Conclusions 


We have presented detailed, multiwavelength, high-cadence 
observations of the unprecedented Type Ic SN 2022jli. We 
attribute the long rise to maximum as the signature of a large 
ejecta mass (Mej ~ 12 + 6 Mo). Future nebular phase spectrosc- 
opy may provide an independent estimate of the core mass from 
the [Ca rr] and [O 1] line ratios (Fransson & Chevalier 1989). 

We provide the first unambiguous detection of periodic 
behavior in an SN optical light curve, measuring a period of 
7.12.5 days and amplitude ~1% of the SN maximum light, 
repeating over a time window of at least ~200 days. This could 


10 


Moore et al. 


be explained by discrete episodes of shock heating from 
interaction with a structured CSM produced through modulated 
mass loss of the progenitor star in a binary system. We also 
consider companion—compact object interaction as the energy 
source but favor a structured CSM. 

We also observe a prolonged early excess and consider two 
scenarios: CSM interaction and ejecta—companion interaction. 
Based on the methods presented in this work, we cannot 
distinguish between these two scenarios. A dense CSM shell 
requires several solar masses of material around the progenitor 
star, requiring exotic phenomena like pulsational pair instabi- 
lity (Woosley 2017) shortly before explosion. Although only 
visible in 10% of cases, we cannot rule out ejecta companion 
interaction, especially given that binarity is already invoked to 
explain the periodic undulations. However, this scenario has 
strict requirements on explosion energy and binary separation. 

SN 2022jli is the subject of further study and multi- 
wavelength observations (Moore et al., in preparation). Late- 
time high-resolution JWST or Hubble Space Telescope 
photometry may reveal the origin of the CSM or a surviving, 
inflated companion star (Liu et al. 2015; Hirai et al. 2018). 
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Facilities: NTT, UH:2.2m, PO:1.2m, Liverpool:2m, LCOGT, 
Gaia, Swift. 

Software: Astropy (Astropy Collaboration et al. 2013, 2018, 
2022), Numpy (Harris et al. 2020), Matlpotlib (Guillochon 
et al. 2018), Mosfit (Guillochon et al. 2018), Hoki (Stevance 
et al. 2020), PSF (Nicholl et al. 2023). 


Appendix 
MOSFiT Posterior 


We the show the posterior distribution in Figure 5 for the 
°°Ni model fit to SN 2022jli described in Section 4.1.2. 
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Figure 5. Posterior distribution for the physical parameters of the MOSFiT default model fit to SN 2022jli. In this model fy; is the fraction of ^5Ni in the ejecta and Maj 
is the mass in solar masses. Ejecta velocity is Vej and log ny,host is the logarithm of the host H column density. Tmin is the temperature floor as defined in Nicholl et al. 
(2017), and t, determines the explosion epoch relative to the earliest photometry point. The c parameter is a white-noise parameter that, when added to all data, gives 


a X? equal to 1. 
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